Momentum distribution of an interacting Bose-condensed gas at finite temperature 
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We use a semiclassical two-fluid model to study the momentum distribution of a Bose-condensed 
' gas with repulsive interactions inside a harmonic trap at finite temperature, with specific focus on 

\ atomic hydrogen. We give particular attention to the average kinetic energy, which is almost entirely 

, ^ . associated with the thermal cloud. A non-linear dependence of the kinetic energy on temperature is 

' displayed, affording a precise way to assess the temperature of the gas. We also show that the kinetic 

energy increases with the strength of the interactions, reflecting an enhanced rate of depletion of 
the condensate with increasing temperature. 



PACS numbers; 03.75.Fi, 05.30.Jp 



I. INTRODUCTION 



The observation in momentum space of Bose-Einstein condensation in a trapped gas with repulsive interactions 
presents at least in principle some advantage over observation in coordinate space. Since the condensate belongs to 
a single quantum state, a broadening of its density profile implies a narrowing of its momentum distribution. At the 
I same time the width of the momentum distribution of the thermal cloud is in essence unchanged, the role of quantal 
. indeterminacy being negligible for this component. 

' The momentum distribution of a dilute Bose-condensed gas has been measured both on trapped and expanding 
. vapours of ^'^Na at T ~ and on trapped vapours of spin-polarized atoms at finite temperature The former 
^ ■ experiment has been successfully interpreted by means of a Thomas- Fermi model |^,^. In the present work we are 
O ' concerned with the latter experiment, in which the momentum distribution of the hydrogen cloud was inferred from 
I the Doppler-sensitive part of the 1S-2S two-photon absorption spectrum: the measured distribution shows a small 

condensate peak on top of a broad thermal component. Cold collisions in the high-density condensate broaden the 
T-H , peak associated with it ||^,^ . Here we are instead interested in extracting thermodynamic information on the system 
^ from the thermal part of the spectrum. 

We use a semiclassical two-fluid model, which was earlier developed to account for the release energy in expanding 
^ j ■ clouds of ^^Rb as a function of temperature [Q. With this model we evaluate the momentum distribution of an in- 
\ teracting Bose gas at equilibrium under harmonic confinement. The results for a choice of parameters corresponding 
■ to the experiment on hydrogen agree with the measured condensate fraction within its error bars. We place partic- 
ular attention on the average kinetic energy of the gas, displaying its temperature dependence and the role of the 
. interactions. 

Before proceeding we should mention that, in addition to the Doppler-sensitive part of present interest, the measured 
spectrum also includes a Doppler-free contribution. The thermal portion of this Doppler-free spectrum has been 
y ' interpreted within an out-of-equilibrium picture invoking large fluctuations in the particle density ^J. However, the 
\ contribution of such fluctuations to the Doppler-sensitive thermal spectrum is expected to be of the same order as for 
^ • the Doppler-free line, making it smaller than the Doppler broadening by two orders of magnitude. 

o 
o 



II. THE MODEL AND ITS JUSTIFICATION 



In our description of the Bose gas with repulsive interactions, confined in a harmonic trap at finite temperature T, we 
^ adopt the semiclassical Hartree-Fock (HF) scheme for the thermal cloud and the Thomas- Fermi (TF) approximation 
for the condensate. 

The choice of a mean field theory is justified by the fact that in current experiments the gas is in a very dilute 
regime, corresponding in the trap to the condition N^^^a/aho 1 @ where N is the number of particles, a is 
the s-wave scattering length and aho = (h/muj)^^^ is the harmonic oscillator length. Here uj — {ui\uiz 

is the 

geometrical average of the harmonic oscillator frequencies characterizing the axially symmetric confining potential 
Vextij) — m{uj\r\ + u)1z^) /2. In particular, in the hydrogen experiment N^/^a/aho — 0.001. 

In the presence of confinement the HF scheme has been shown to provide a simple and accurate description of the 
excitation spectrum []l0|], yielding results which are comparable with the predictions of the full Bogoliubov theory. 
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The semiclassical approximation is valid for a bosonic cloud when its temperature is appreciably larger than the 
harmonic-oscillator spacing, fc^T ^ huj±^z- In particular, in the hydrogen experiment ksT ~ 270huj± ~ lO^hcUz. For 
the specific case of trapped atomic vapours, the semiclassical HF approximation has also been tested by a Monte 
Carlo calculation performed with a choice of parameters corresponding to the ^^Rb experiments |^ . 

The TF approximation for the condensate well describes the gas with repulsive interactions in the strong coupling 
limit Nca/cLho S> 1, *-e. when the number Nc of atoms in the condensate is large. This approximation is therefore 
valid everywhere except in the critical region. In the hydrogen experiment one has Nca/aho — 1-6 • 10'*. 

On all these counts, therefore, the use of the semiclassical two fluid model is fully justified in the calculations 
that we present below, except in the immediate neighbourhood of the critical temperature as remarked just above. 
The momentum distributions of the thermal cloud /t(p) and of the condensate /c(p) are then given by the Bose 
distribution in an effective potential and by the square modulus of the Fourier transform of the wavefunction of the 
condensate: 
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Here </>c(r) ~ \/nc{r), the square root of the density profile of the condensate. In Eq. (||) the effective potential acting 
on the thermal cloud is given by 



Veff{r) = Vea:t{r) + 2gncir) + 2gnT(r) , 
nT(r) being the density profile of the thermal cloud. Finally, we have 

nc(r) = i - Vextir) - 2gnT{r)] e{fj. - Vext{r) - 2gnT{r)) 

and 
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The coupling constant g is fixed by the s-wave scattering length a through the relation g — Anti^a/m, and the chemical 
potential /i is fixed by the total number N of atoms, 



J [n,{r)+nTir)]dh = N 



(6) 



For temperatures above the critical condensation temperature the expression (|l|) is equivalent to that developed by 
Chou et al. ill. 



III. MOMENTUM DISTRIBUTIONS 

In the experiments of Fried et al. j^] on atomic hydrogen the momentum distribution has been inferred from the 
Doppler-sensitive spectrum as measured along the weakly confining axis of the trap. We accordingly report in Figures 
1 and 2 the distributions for the thermal cloud and for the condensate as functions of the axial momentum pz : 

FtAp^) - / ^2fTAv)- (7) 

The units of length and momentum are ahoL = {h/'mijj±_Y/'^ and phoi. = {limuj±Y^'^ , with the values of m and uj±_ of 
the experiment [Q. We have taken uj± = 2tt ■ 3.9 kHz and = 27r • 10.2 Hz. 

Figure |l| illustrates both the role of statistics and the effect of the interactions on the axial momentum distribution 
of the thermal component. The predictions of the semiclassical model for a weakly interacting gas, with a scattering 
length chosen equal to that of ^H {qh = 6.48 • 10^^ nm from accurate calculations by Jamieson et al. [|l3| ) or arbitrarily 
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increased to 20 au for illustrative purposes, are compared with each other and with that for a non-interacting Bose 
gas. The value of the temperature is T = 50.3 /LtK, as will be discussed in Sect. below. Figure |l| also reports the 
momentum distribution of the classical gas at the same temperature. 

It is evident from Figure ^ that quantum statistics increases the population of the states at low momentum, as 
expected. The main effect of the interactions is instead a depletion of the condensate, leading to increases in the 
population and in the width of the momentum distribution of the thermal cloud. 

The momentum distribution of the condensate at the same temperature as coming from the Thomas-Fcrmi model 
is shown in Fig. |^. 

IV. KINETIC ENERGY 

We proceed to study the momentum distribution of an interacting Bose gas as a function of temperature by focusing 
on its second moment, that is on the kinetic energy of the trapped cloud. This is the quantity which can be most 
directly compared with the experimental data . 

In Figure ^ we show the behaviour of the kinetic energy of the gas for two values of the scattering length, correspond- 
ing to the experimental value a = an for the hydrogen gas and to the stronger interaction a = 20aH, as compared 
with the non-interacting quantum and classical results. For the parameters corresponding to the hydrogen experiment 
the effect of the interactions on the kinetic energy is clearly quite negligible. Indeed, the dilution parameter for this 
system is very small, N^^^an / aho — 0.001 as already noted. On increasing the interaction strength we find that the 
kinetic energy increases. This effect can be understood as coming from enhanced depiction of the condensate. There 
is instead a major effect of quantal statistics on the average kinetic energy below the condensation temperature, as is 
seen from Figure |[ 

It is important to stress that in the interacting gas the behaviour of the kinetic energy depends not only on the 
scattering length but also on the total number of particles in the system. This is illustrated in Figure ^. We have 
exploited this fact to test how far the semiclassical model can quantitatively explain the experimental results of Fried 
et al. on the hydrogen gas. After choosing a value for the total number of particles in the cloud within the 
experimental error bars, we have fitted to experiment the calculated kinetic energy and extracted from the model the 
temperature of the cloud and the condensate fraction. For a total number of atoms equal to TV = 10.5 • 10^ we obtain 
Nc/N = 10%, that is = 1.05 ■ 10^ at a temperature T = 50.3 /iK. The experimental values are Nc/N = 6^3% and 
Nc = 1.1 ± 0.6 • 10^. Therefore, we conclude that the semiclassical model is compatible with the experimental data 
for the momentum distribution of the thermal cloud. 

Our calculations also show that from kinetic energy data it is possible, through the semiclassical model to obtain 
a good estimate for the temperature of the Bose cloud below the critical temperature. In the intermediate- to low- 
temperature region this procedure yields a result which differs significantly from the classical result E^^^ = ksT /2. 
For the specific case of the hydrogen experiment we estimate that the value of the temperature corresponding to the 
measured kinetic energy is T = 50.3 /iK, that is about 20% higher than the value extrapolated from the classical result 
(see again Figure ^ . 

V. SUMMARY AND CONCLUSIONS 

In summary, the Doppler-sensitive part of the 15-25 two-photon absorption spectrum of Bose-condensed atomic 
hydrogen contains quantitative information on thermodynamic properties, which can be extracted from the available 
data by means of an analysis combining a Thomas-Fermi treatment of the condensate with a semiclassical Hartree- 
Fock treatment of the thermal cloud. The spectral broadening due to fluctuations in the density is negligible and the 
spectral width is directly related to the kinetic energy of the gas. 

Within such a semiclassical two-fluid model we have studied the momentum distribution and the mean kinetic 
energy of a Bose-condensed gas with repulsive interactions. From the behaviour of the kinetic energy as a function 
of temperature and of the coupling strength we have reached two main conclusions: (i) the interactions enhance the 
rate of depletion of the condensate with increasing temperature, causing an increase of the kinetic energy; and (m) 
the relation between kinetic energy and temperature is strongly non-linear, owing to the role of quantal statistics as 
opposed to classical statistics. It appears that the present model may yield useful improvements in the analysis of 
experimental data. 
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FIG. 1. Axial momentum distribution of the thermal component of a weakly interacting Bose gas with N = 10.5 • 10^ at 
T = 50.3 /jK. The scattering length is a = 20aH (dotted curve) and a = an (continuous curve) with uh = 6.48 • 10~^ nm. These 

two curves arc compared with the momentum distribution of the non-interacting Bosc gas (dashed curve) and of the classical 
gas (long-dashed curve). The latter has the same number of particles as the non-interacting quantal gas. 
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FIG. 2. Axial momentum distribution of the condensate in an interacting Bose gas with N = 10.5 • 10^ at T = 50.3/xK. The 
scattering length is a = an- 
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FIG. 3. Kinetic energy versus temperature (in units of the ideal-gas critical temperature To) for an interacting Bose gas with 
N — 10.5 ■ 10^ at a scattering length a — 20a_r/ (stars) and at the hydrogen scattering length (filled circles). These results are 
compared with the kinetic energy of the non-interacting Bose gas (dashed curve) and of the classical gas (dotted curve). The 
horizontal full line locates the kinetic energy measured in the hydrogen experiment M. 
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FIG. 4. Kinetic energy versus temperature for an interacting Bose gas at different values of the total number of particles: 
N = 18.3-10^ (long-dashed), N = 15-10^ (dashed-dotted), N = 11.6-10^ (dotted), iV = 10.5-10*' (continuous) and N = 9.5-10^ 
(short-dashed). The horizontal full line locates the kinetic energy measured in the hydrogen experiment Q|. 
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